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Lensing 
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Abstract. The lensing information provided by multiple images, ar- 
clets, and statistical distortions can all be formulated as linear constraints 
on the arrival-time surface, and hence on the mass distribution. This re- 
\ duces cluster lens reconstruction from combined strong and weak lensing 

0^ ■ to a standard type of inversion problem. Moreover, the mass sheet de- 

CD 1 generacy is broken if there are sources at different redshifts. 

This paper presents a reconstruction technique based on these ideas, 
and a mass map and discussion of CI 1358+62. 



1. Introduction 



The different regimes of cluster lensing — the inner region where multiple images, 
including giant arcs are found, further out where highly sheared but singly im- 
^ . aged arclets are found, and the outer regions showing statistical shear — have in 

5h I the past been approached with quite different modeling methods. But these ap- 

parently separate regimes can be studied in a unified way. The key is to express 
lensing information in terms of the arrival surface. 

In earlier work (see AbdelSalam et al. 1998) we combined the multiple- 
image and arclet regimes and argued that extending to statistical shear was a 
straightforward algorithmic issue. This paper makes that extension. 

2. The arrival-time surface 

The creature we will mostly be concerned with is the scaled arrival-time surface 
t(9), and it is expressed as follows: 

r (0) = i(e _ _ (As/ y ln |g _ K{Ql) d 2 6 , (1) 

= \(6-(3f -2{D 1s /D s )V- 2 k(G). (2) 
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Figure 1. Schematic of an arclet elongated with position angle u>, 
with (9 x i,8 y i) being a coordinate system aligned with the elongation. 



As usual, 8 and /3 are angular locations on the image and source planes respec- 
tively, and the D's are angular diameter distances in units of c/Hq. However, k 
is the convergence for sources at infinity, (hence the factors of D\ s /D s ) to make 
dealing with multiple source redshifts more convenient. The form of Eq. (|l]) is 
then familiar, while Eq. (||) is just a shorthand for the same thing with V -2 
denoting an 'inverse Laplacian operator' in two dimensions. 

To get back to physical arrival time and surface density we use 

time = h~ x z x 80 days arcsec -2 x t(9) (3) 
surface density = h~ l z x 1.2 x lO n M arcsec -2 x k(0) (4) 

where z = (1 + z\)D\ but is ~ z\. 

Lensing data constrain the arrival-time surface and hence the mass distri- 
bution k(8) in various ways. We can identify three types of constraints: on the 
values of the arrival time surface at particular points, on the first derivative, and 
on the second derivative. 

The first type is when we know the height difference between two or more 
points on the arrival-time surface. This is all-important in time-delay quasars 
and Hq measurements, because time delays supply just this height-difference 
information. But it is not relevant to cluster lensing, at least not yet, though 
we can always hope for a supernova in a giant arc. 

The second type of constraint comes from stating Fermat's principle at each 
image location, i.e., 

Vr(0 image ) = 0. (5) 

Since image positions ©image can generally be measured very accurately, the 
Eq. (|B|) tells us that the arrival-time surface has zero gradient at some known 
point. This type of constraint is only useful if we have multiple images; in that 
case we have 2 (images) equation of the form (||) but only two unknown source 
coordinates to solve for, giving us some net constraints on k{8). In general there 
are 2((images) — (sources)) constraints. 

The third type of constraint involves the curvature or second derivative of 
the arrival-time surface, and may come from either arclets or from statistical 
distortions. 

Consider first a situation like Fig. 1, where an arclet is observed elongated 
with position angle u. If we are confident that the arclet has been elongated by 
a factor of at least k, we can consider the coordinate system (8 x /,8 y i) aligned 
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with the elongation and write 



d 2 

d9 2 , T ^ image ' 
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Qq2 M^image, 



(6) 



Further, since u> is known, we can transform Eq. (0) to the unrotated (9 x ,9 y ) 
coordinates. It is not necessary to have highly accurate values of the elongation 
or its orientation oj — all we have to do is set k conservatively enough that the 
inequality (||) is valid. The arclet used may itself be part of a multiple-image 
system; if so we will have to guess the parity to remove the absolute value signs 
in®. 

The same idea can be applied to statistical distortion. In this case the data 
provide an estimate and uncertainty 
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9i ± Aft, 
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where ji and k are components of the second derivative of t{6). We can rewrite 
this in various ways, such as 

(1 - K(0image))9i ~ < Tutelage) < (1 ~ /«(^imago))5j + ^9i- (8) 

(We have assumed nAgi is negligible here.) 

A key feature of the constraint equations (||), (|6|),.-and (||) is that they are 
all linear in the unknowns (3 and n{0). Which is to sayu 



Lensing 
data 
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(9) 



and moreover this is true for all lensing regimes: multiple-image, arclet, and 
statistical shear. 



3. Mass reconstruction 

Eq. (^) summarizes the advantage of casting the observational information as 
constraints on the arrival-time surface. Lens reconstruction is reduced to a linear 
inversion problem. It is (as symbolized by the 2x5 matrix) a highly underde- 
termined problem, so in order to reconstruct the lens it is necessary to add extra 
information (sometimes called a prior) . Also as indicated in the 2x5 matrix, the 
data enter into the linear operator — an unusual complication. However, as we 



1 With one exception: scalar magnification data are quadratic in k and so (g) does not apply. 
Taylor and collaborators (see e.g., Dye & Taylor 1998) have developed reconstruction techniques 
for this case. However, if magnification information is present along with shear information, 
together they give (at least in principle) the tensor magnification which is linear in k, and so 
(IpJ) again applies. 
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Figure 2. At fixed redshift, the mass sheet degeneracy leaves 1 — k 
uncertain by a multiplicative factor, while the monopole part of k is 
undetermined outside the region of the images. As a result, the two 
mass profiles shown will produce identical multiple-image and shear 
data out to the cliff. 



discussed above, the data enter there either as image positions which are known 
very accurately, or as inequalities which can be set conservatively; so this issue 
does not introduce new difficulties. 

At least four different possibilities for lens reconstruction now suggest them- 
selves 

• Put r{0) on a grid and regularize. This may be the simplest approach. 

• Use basis function expansions for k{9) and t(0). A Fourier-Bessel expan- 
sion 

T ( e ) = ^2 c mnJm(k mn 9) exp(im(f)) 

mn 

is particularly attractive as it would trivially relate r and K. 

• Pixellate k(0) and regularize. This is the one we have implemented. The 
mass is distributed on square tiles, each having constant but adjustable 
density, and t(0) comes from computing Eq. (0) exactly. Note that al- 
though the mass distribution is discontinuous, the arrival-time surface is 
smooth. Unlike in the two previous possible approaches, pixellating the 
mass makes it is easy to enforce non-negativity of the mass distribution. 

• Pixellate k(0) and use maximum entropy. See Bridle et al. (this volume) 
for an example of this in a somewhat different context. 

The regularization we applied is to minimize 

J (k- 2 d 2 + e 4 J (V 2 k) 2 d 2 (10) 



while of course enforcing all the lensing constraints. The first term in (10) tends 
to minimize mass-to-light variation, since (k)/{L) is a mean M/L. The second 
term tends to smooth, with e a sort of smoothing scale. By regularizing with 
respect to different light distributions we can get an estimate of the uncertainty. 
Since the regularizing functional (|i~0|) is quadratic in k, minimizing it sub- 



ject to the linear lensing constraints is readily implemented through quadratic 
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Figure 3. Left: Cluster luminosity contours in CI 1358+62. from 
the catalog of galaxy VzS in van Dokkum et al (1998). Likened at 
the conference to a left facing seahorse with a ponytail. Right: Shear 
map (exaggerated by a factor of 5) from Hoekstra et al (1998). We 
have excised the central region to avoid averaging the shear where it 
is varying rapidly. Also note that unlike Fig. 14 in Hoekstra et al., we 
have binned without smoothing here. 

programming algorithms. The disadvantage is that the storage needed is ~ 
2 (number of pixels) 2 limiting us to of order 5000 pixels. The solution is to have 
adaptive pixellation (much like tree codes in iV-body simulations), with smaller 
pixels for the inner parts of the cluster and large pixels outside. 




4. The mass-sheet degeneracy 

After enthusing about how easy lens reconstruction is, it is well to make a 
cautionary remark about the main source of uncertainty. Here again, the arrival- 
time surface is very useful. 

If we multiply r(0) by a constant factor, we just stretch the surface verti- 
cally, and geometrically it is clear that neither image locations nor their relative 
magnifications change. More formally, we can rearrange Eq. (0) and drop the 
irrelevant term ^/3 2 to get 

r(0)=2V- 2 (l-^)-±0-/3. (11) 

If we now multiply both (1 — D\ s /D s k) and (3 by some constant a, the image 
structure will be unchanged; only time delays and total magnification will get 
multiplied by a and a~ l respectively (the latter because rescaling f3 rescales 
all the sources). Therefore lens reconstruction from image structure (without 
absolute magnifications) leaves (1 — D\ s /D s k) uncertain by a constant factor. 
This is the mass-sheet degeneracy. 
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Figure 4. Pixellation used for the mass reconstruction of the Sea- 
horse cluster. Note that is completely adaptive, and allows us to 
increase the resolution at any desired place — in this case the region 
around the red arc. Left: The full field. Right: The inner region, 
with dots marking features in the red arc. 



It may appear at first that the mass-sheet degeneracy can be eliminated by 
requiring n to vanish at large distances. But outside the observed region, the 
monopole part of k is completely unconstrained. So the mass-sheet degeneracy 
is equally effective with a mass disk larger than the observed region. Figure 2 
illustrates. 



5. Reconstruction of CI 1358+62 (The Seahorse Cluster) 

This cluster seemed the obvious test case for combined strong and weak lensing 
reconstruction in a relatively large field. The cluster itself is at z = 0.33 and its 
inner part lenses a z = 4.92 galaxy into a red arc; this was identified by Franx 
et al. (1997) who also presented a strong lens model. Hoekstra et al. (1998) 
measured the shear field in an HST WFPC2 mosaic and reconstructed a mass 
map from weak lensing. The present work combines both regimes. 

Figure 3 illustrates the data on this cluster that we have used: a shear field 
derived by binning the individual polarization measurements kindly provided by 
Hoekstra and collaborators and assuming a constant z = 1 for the background 
galaxies, and a smoothed V z light distribution. Figure 4 shows the pixellation 
we used. 

Figure 5 is our mass map, computed by regularizing with respect to the 
light and a smoothing scale e changing from ~ 5" in the center to ~ 1' at the 
edge. Our estimated uncertainty is derived from an ensemble of reconstructions 
where we rotated the light map by arbitrary angles and shifted it randomly by 
up to 100" and regularized with respect to these altered light maps. The mass 
map resembles Fig. 15 of Hoekstra et al. (1998) but tends to be smoother. Also 
the overall normalization is somewhat higher; this is probably due to different 



6 




Figure 5. Mass map of the Seahorse cluster. Left: Contours of k in 
steps of 0.1. Right: Contours of An in steps of 0.025. 



treatments of the boundary, though with the mass-sheet degeneracy. (The red 
arc at z = 4.92 compared to z = 1 assumed for the weakly lensed galaxies, 
reduces the effect of the mass-sheet degeneracy, but does not eliminate it.) Our 
central density is higher, but this is as expected since inclusion of a multiply- 
imaged system immediately forces k > 1. There is some indication that the 
mass peak is offset (by some 10s of kpc) to the south of the light peak, a natural 
thing to expect if the cluster is asymmetric and galaxy formation is biased. But 
this offset is tentative; in Abell 2218, where the inner region is much richer in 
lensing and better constrained, the evidence for an offset is more compelling 
(AbdelSalam et al. 1998). 

Figure 6 shows the enclosed mass out to different radii and the reconstructed 
shear. It is surprising that the enclosed mass after angular averaging, looks so 
'isothermal' even though the mass map is very asymmetric. We obtain M = 
(10 ± 1) x lO 14 M0Mpc _1 , corresponding to a formal Einstein radius of 43" 
and a formal isothermal los dispersion of (990 ± 70) km sec -1 . The estimated 
mass-to-light is (380 ± 60) /i M L~^. The CNOC survey (Carlberg et al. 1997) 
measured 910 km sec -1 for the los velocity dispersion and estimated M/L = 
229h- 1 M P) L-} / . 



6. Discussion 



We have developed a mass reconstruction technique that deals with all the clus- 
ter lensing regimes simultaneously, from multiple images in the central regions 
to weak shear in the outer regions, and moreover with adaptive resolution. Code 
implementing this is available from the authors. Variants of our technique (sug- 
gested in Section 3) may also be of interest in future work. 

At this stage it appears that cluster mass reconstruction is very good at 
recovering features — a good example are offsets between mass and light peaks 
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Figure 6. Left: Shear map. The superimposed ellipses show the ob- 
served and reconstructed (dotted) shear, both exaggerated by a factor 
of five. As we might expect, the reconstructed shear appears more cor- 
related than the observed. The ellipses in the inner region show the 
reconstructed shear (not exaggerated) at the locations of the red arc. 
Right: Enclosed mass and estimated uncertainty. 



in clusters, which may be indirect evidence for biased galaxy formation. But 
the calibration of the mass (even with perfect shear data) is still problematic — 
the mass sheet degeneracy is a vicious effect. We suggest that recovering the 
redshift distribution and the absolute magnification of the background galaxies 
(even with large uncertainties) are the best hope of breaking this degeneracy, 
and is likely to be very rewarding in future work. 

Acknowledgments. We are grateful to Henk Hoekstra and collaborators 
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